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ABSTRACT 

We performed high-resolution simulations of two stellar collisions relevant for 
stars in globular clusters. We considered one head-on collision and one off-axis collision 
between two 0.6 Mq main sequence stars. We show that a resolution of about 100 000 
particles is sufficient for most studies of the structure and evolution of blue stragglers. 
We demonstrate conclusively that collision products between main-sequence stars in 
globular clusters do not have surface convection zones larger than 0.004 Mq after 
the collision, nor do they develop convection zones during the ‘pre-main-sequence’ 
thermal relaxation phase of their post-collision evolution. Therefore, any mechanism 
which requires a surface convection zone (i.e. chemical mixing or angular momentum 
loss via a magnetic wind) cannot operate in these stars. We show that no disk of 
material surrounding the collision product is produced in off-axis collisions. The lack 
of both a convection zone and a disk proves a continuing problem for the angular 
momentum evolution of blue stragglers in globular clusters. 

Key words: blue stragglers - globular clusters - hydrodynamics - convection - stars: 
rotation 


1 INTRODUCTION AND MOTIVATION 


It has been shown that the products of main-sequence - 
main-sequence collisions appear in the colour-magnitude 
diagrams of clusters a s blue stragglers ( Sills et al. 1997 
Sandquist et al. 1997). Since blue stragglers are readily 
observable in clusters, they form an ideal population with 
which to probe the dynamical evolution of the cluster. The 
dynamical state of a globular cluster (its density profile, ve¬ 
locity dispersion, amount of mass segregation etc.) will de¬ 
termine the rate and nature of the collisions which occur 
in the cluster. As the cluster evolves, the kinds of collision 
that occur will change. Therefore, the population of colli¬ 
sion products in a cluster can be used to probe the history 
of the cluster (Sills et al. 200C; Hurley et al. 2001). However, 
in order to use collision products in this way, there are two 
issues which must first be understood. First, we know that 
blue stragglers can also be formed through the merger of two 
components of a binary system. These blue stragglers will 
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probably have different properties than those formed from 
collisions. We must either be confident that the population 
we are observing is collisional in origin (e.g. from cluster 
density considerations), or be able to distinguish between 
the two populations. Secondly, we also need to be sure that 
we understand the formation and evolution of the collision 
products themselves. This paper is concerned with address¬ 
ing the second point. 


In this paper, we present the highest resolution 
smoothed particle hydrodynamic (SPH) simulations of col¬ 
lisions between main-sequence stars to date. Most recent 
computations have ~ 10^ particles (Lombardi et al. 1996; 
Sandquist et al. 1997), wi th the highest re solution simula- 
tion using 10’” particles (sills et al. 2001). In this paper, 
we increase the number of particles to 10“. There are three 
main reasons to extend this kind of simulation to such high 
resolution. The first, and simplest, is to make sure that no 
fundamental change in our understandi ng of blue stragglers 
occurs. In earlier work, Benz & Hills (1987) performed an 
SPH simulation using 1024 particles, and concluded that 
collision products are fully mixed (i.e. the resulting star 
is chemically homogeneous). Ten years later, these simula¬ 
tions were repeated but with a factor of 10-50 more particles 
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(Lombardi et al. 1996). Because of the higher resolution, it 


became clear that collision products are NOT chemically ho¬ 
mogeneous, but instead retain some memory of the chemical 
profile of the parent stars of the collision. We feel conhdent 
that such a fundamental change in our understanding will 
not happen again when we increase the resolution, but we 
want to be sure. 

More importantly, we performed these high resolution 
calcula tions to answer two fundamental questions about the 


structu ’e of the collision products immediately after the col- 
lision. The first involves the structure of the outer layers of 
the star, and the second is concerned with the material which 
is thrown off by the star during the collision itself. There has 
been some debate in the literature recently about whether 
the collision product has a surface convection zone shortly 
after the end of the collision. The presence of such a con¬ 
vection zone could have ramifications for both the surface 
chemical abundances and the rotation rate of the star when 


We use the standard form of artificial viscosity with a = 1 
and (3 = 2.5, and an adiabatic equation of state. The ther¬ 
modynamic quantities are evolved by following the change in 
internal energy. Both the smoothing length and the number 
of neighbours can change in time and space. The smoothing 
length is varied to keep the number of neighbours approxi¬ 
mately constant (~ 50 for the low resolution runs, and ~ 200 
for the million particle run). The code is the parallel version 
of th e code described in detail in Bate, Bonnell & Price ( 
1995). It was parallelized using OpenMP, and was run on the 


it reaches its main sequence (Leonard & Livio 1995). Previ¬ 
ous SPH simulations and subsequent evolution calculations 
have shown that no convection zone exists at the end of the 


SGI Origin 3000 operated by the UK Astrophysical Fluids 
Facility (UKAFF) based at the University of Leicester. 

We modelled collisions between two equal-mass stars. 
We chose a mass of 0.6 Mq as a representative, but not 
extreme, mass for globular cluster stars. The initial stel¬ 
lar models were calculated using the Yale stellar evolution 
code (YREC, Guenther et al. 1992), and had a metallicity 
of Z=0.001 and an age of 15 Gyr. If the goal of the simula¬ 
tion is to investigate the detailed structure of the collision 
product, then it is crucial to begin with a realistic st ellar 
model rather than a polytrope or other approximation (Bills 


collisioif, nor does it appear during the initial thermal re- _ fc Lombardi 1997|). The particles were initially distributed 


taxation of the star (Lombardi et al. 1996; Sills et al. 1997). 
However, members of the community argue that the previ¬ 
ous simulations are unable to resolve the outermost regions 
of the star well at all, due to their low particle number and 
use of equal-mass particles (Livio, private communication). 
In this paper, we resolve this issue by increasing the reso¬ 
lution of our simulations in the outer regions of the parent 
stars in particular, and by using variable mass particles. 

We are also interested in following the material which 
is thrown off from the stars during the collision. We have 
discovered, using previous generation simulations, that blue 
stragglers which are formed by an off-axis coll ision have 
an angular momentum problem ( ^ills et al. 2001 ). Namely, 
these stars retain too much of their angular momentum, and 


on an equally-spaced grid, and their masses were varied un¬ 
til the density profile matched that of the stellar model. By 
using unequal mass particles, we increase the resolution in 
the outer, low density regions of the star. This is important 
for resolving the outer layers of the collision products (to de¬ 
termine if there is a convection zone or not) and for following 
the material which is thrown off during the collision. 

The stars were given a relative velocity at infinity of 
Voo = 10 km/s, which is a reasonable value for globular cluster 
stars. They were set up on almost parabolic orbits with a 
pericentre separation rp = 0 for the head-on collisions, and 
rp=0.25 i ?0 for the off-axis collisions. (The radius of the 
main sequence star is 0.617?©.) Most collisions are ex pected 


have nC' apparent way of losing it during their thermal re¬ 


laxation phase after the collision. As a result, the collision 
products spin up during their collapse to the main sequence, 
and inevitably rotate faster than their break-up velocity. 
Such collision products can never become blue stragglers, 
since they tear themselves apart before they reach the main 
sequence. We wish to follow the lost material with greater 
resolution, to study the amount of angular momentum that 
this material carries away with it. We also wish to follow 
the outermost material of the parent stars to see if some of 
it forms a disk around the collision product. If so, then we 
can plausibly suggest that magnetic locking to a disk is re¬ 
sponsible for removing angular momentum from the collision 
product. 

In ^ we describe the methods used to model the stellar 
collisions and their subsequent evolution. We present our 
results in and discuss their implications in §^. 


2 METHOD 

The simulations of stellar collisions discussed in this paper 
were performed us ing the smoothed particle hy drodynam- 
ics (SPH) method (Benz 199C; Monaghan 1992). Our three 
dimensional code uses a tree to solve for the gravitationa l 
forces and to find the nearest neighbours (Benz et al. 1990). 


pec 

to happen du ring interactions involving binary stars (|Hut & 


Bahcall 1983), so the stars are on bound, and usually highly 


elliptical, orbits. For both the head-on and off-axis collisions, 
the number of particles was varied between 10 000 and 300 
000. We also ran one head-on collision with 10® particles. 
The details of the collisions are shown in Table |^. We give 
a single letter name for each run, the pericentre separation 
rp in solar radii, the number of particles, the initial value 
of the smoothing length ho in solar radii, the minimum and 
jmaximum particle masses in solar masses, and the amount 
of mass lost from the system during the collision in solar 
masses. 

The results of these three dimensional SPH simulations 
were converted to one dimensional models suitable for start¬ 
ing models for the Yale stellar evolution code, YREG. The 
entropy of the particles and their chemical abundances were 
averaged over surfaces of constant gravitational potential 
and binned into ~ 100 bins. The structure of the collision 
product is determined from these profiles using the equation 
of hydrostatic equilibrium. The temperature profile is calcu¬ 
lated using the ideal gas equation of state. For more details 
of the conversion between SPH results and stellar evolution 


models, see Sills et al. (1997; 2001). 

We used YREC to calculate the evolution of the colli¬ 
sion products. Using the starting models described above, 
we followed their evolution during the thermal relaxation 
phase (analogous to the pre-main-sequence of normal stars). 
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Table 1. Description of SPH simulations performed 


Run 

Tp 

Rq 

-^part 

ho 

Rq 

M© 

M© 

Most 

M© 

A 

0.0 

10 162 

0.058 

2.10 X 10 “® 

4.58 X 10 -® 

8.14 X 10 -° 

B 

0.0 

29 950 

0.0417 

3.71 X 10 -'^ 

1.37 X 10 -® 

7.76 X 10 -° 

C 

0.0 

100 034 

0.028568 

4.96 X 10 “® 

4.47 X 10-1 

7.49 X 10 -° 

D 

0.0 

299 398 

0.02011 

5.27 X lO ”'^ 

1.56 X 10-1 

7.41 X 10 -° 

E 

0.0 

999 778 

0.0136 

2.98 X 10 - 1 ° 

4.75 X 10 -® 

7.47 X 10 -° 

F 

0.25 

10 162 

0.058 

2.10 X 10 “® 

4.58 X 10 -® 

3.45 X 10 -° 

G 

0.25 

29 950 

0.0417 

3.71 X 10-1 

1.37 X 10 -® 

3.62 X 10 -° 

H 

0.25 

100 034 

0.028568 

4.96 X 10 -® 

4.47 X 10-1 

3.99 X 10 -° 

I 

0.25 

299 398 

0.02011 

5.27 X 10 -° 

1.56 X 10-1 

4.16 X 10 -° 


through their main-sequence lifetimes and to the giant 
branch. 


3 RESULTS 

We performed five simulations of head-on stellar collisions, 
with the number of particles ranging from lO'* to 10®. The 
average resolution elements for these collisions ranged from 
2h = O.1167?0 to 2h = 0.0277?© respectively. Figure ^ shows 
the profiles of the important structural parameters as a func¬ 
tion of mass fraction in the product of a head-on collision 
between two 0.6 M© stars. The different line styles denote 
increasing particle numbers in the following order: 10 000 
(solid), 30 000 (dotted), 100 000 (short dashed), 300 000 
(long dashed), 1 000 000 (dot-dashed). The low resolution 
simulations do a reasonable job of producing the correct 
structure for the collision products. However, as the reso¬ 
lution increases, the models converge to an answer differ¬ 
ent from that seen in the low resolution runs. This is best 
seen in the plot of hydrogen abundance as a function of 
mass fraction. The 100 000 and 300 000 particle simulations 
have converged to a solution for the structure of this colli¬ 
sion product. The higher resolution simulations also produce 
smoother profiles. This is expected, since we have the infor¬ 
mation from more particles contributing to the values of the 
structural quantities (pressure, density, etc.) in each shell. 
Therefore, we suggest that the optimum number of parti¬ 
cles for studying the structure of stellar collision products is 
around 100 000. 

In Figure ^ we present the evolutionary tracks which 
result from the starting models shown in Figure The stars 
begin their post-collision lives at position A. They are quite 
bright, and are significantly out of thermal equilibrium be¬ 
cause of the energy deposited into their envelopes by the 
collision. They contract to the main sequence (position B), 
where core hydrogen burning begins. The stars look like nor¬ 
mal main-sequence stars, except that their central hydrogen 
abundance is less than it would be for a normal star of that 
mass. Therefore, their main-sequence lifetime is shorter than 
that of a normal 1.1 M© star (3.3 Gyr). The collision prod¬ 
ucts reach the turnoff, where hydrogen is exhausted in their 
core, and evolve towards the giant branch (position C), ex¬ 
actly like normal stars. The time since the collision at three 
benchmark positions along the evolutionary track is given 
in Table |^. Both the post-collision structure and the evolu- 






Figure 1. Structure of head-on collision products. The main 
structure parameters of a non-rotating star have been plotted as 
a function of mass fraction: pressure, temperature, density, and 
hydrogen mass fraction. The different line styles show the results 
for different numbers of SPH particles. In increasing order of par¬ 
ticles, they are: solid (10 000 particles), dotted (30 000 particles), 
short-dashed (100 000 particles) long-dashed (300 000 particles), 
and dot-dashed (1 000 000 particles). The models converge to¬ 
wards the 1 000 000 particle simulation, and there is very little 
obvious difference between the 100 000 , 300 000 and 1 000 000 
particle results. 


tionary t racks agree with previous simulations of the same 
collision (Sills et al. 1997). 

The different line styles in Figure indicate the dif¬ 
ferent number of particles in the simulation, with the line 
styles having the same meaning as in Figure |^. We see the 
same trend with increasing resolution of the simulation - 
the tracks converge to a common solution around 100 000 
particles. This is expected, since the evolution of a star is 
determined by its structure. The same trend of convergence 
is seen in the evolutionary timescales in Table Note that 
the largest differences between simulations show up in the 
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Table 2. Evolutionary Timescales (Gyr) 


Run 

Return to MS 

Turnoff 

Giant Branch 

A 

1.60 xlO-3 

1.41 

2.08 

B 

1.90 xlO-3 

1.64 

2.41 

C 

2.20 xlO-3 

1.85 

2.69 

D 

2.16 xlO-3 

1.92 

2.79 

E 

2.20 xlO-3 

1.87 

2.73 



4 3.9 3.8 3.7 3.6 


log (T,„) 

Figure 2. Evolutionary tracks for head-on collision products. The 
stars begin their post-collision evolution far from the Hayashi 
track, but are still bright (position A). As they contract, they 
move towards the main sequence (position B), and then follow a 
standard evolutionary track towards the giant branch (position 
C). The line styles are the same as in Figure The difference 
in structure between the different resolution simulations mani¬ 
fests itself as a difference in evolutionary track shape. The tracks 
are most different during their initial contraction, and then again 
near the turnoff of the main sequence. The lowest resolution sim¬ 
ulation (10 000 particles, solid line) has a hook in its evolutionary 
track near the turnoff, evidence of a well-developed core convec¬ 
tion zone. The higher resolution simulations do not show this 
feature. 

thermal relaxation phase, and that by the main sequence, 
all the simulations are almost indistinguishable. The excep¬ 
tion is the lowest resolution simulation. The solid black line 
begins its main sequence quite close to the dotted line (30 
000 particle simulation), but has a hook near the turnoff, 
indicating the presence of a central convection zone. As the 
resolution increases, we see that this convection zone disap¬ 
pears. Simulations with 10^ particles are not quite detailed 
enough to accurately depict the true structure of these col¬ 
lision products. 

The criterion for convective stability in stellar models 
is the Schwarzschild criterion: 

V-Vad<0 (1) 



0 0.2 0.4 0.6 0.8 1 

M/M^ 


Figure 3. Convective stability criterion V — V^d as a function of 
mass fraction for the head-on collision models immediately after 
the collision. The line styles are the same as in Figure ^ When 
V — Vad is less than zero, the region is stable against convection. 
Except for a small region at the centre of the star, the entire 
collision product is stable against convection. 

where Vad is the adiabatic temperature gradient in the star, 
and V is the radiative temperature gradient: 

^ dlnT 3 kLP 
^ din GA-KaG MT* 

where a is the Stefan Boltzmann constant, G is the 
gravitational constant, and k,L,P,M and T are the opac¬ 
ity, luminosity, pressure, enclosed mass and temperature at 
that position in the star. The models taken from the SPH 
results were used as starting models in the stellar evolution 
code, and allowed to relax, so that all the equations of stellar 
structure are satisfied. Since our SPH code does not allow 
for energy transport, the luminosity distribution in the star 
can only be determined from the other structure parameters 
using the equations of stellar structure. The temperature 
gradients were calculated, and their difference is plotted as 
a function of mass fraction in Figure The line styles are 
the same as in Figure We can clearly see a convective 
core in these stars, where V — Vad > 0. However, the entire 
star outside of the inner 0.1 Mq is not convective. The out¬ 
ermost shell in all our stellar models is at a mass fraction 
M/M* = 0.99667. This value is the same for all models, re¬ 
gardless of the number of SPH particles, because of the way 
we have transformed the SPH information to the stellar evo¬ 
lution starting model. The total mass of these stars is 1.12 
Mq, so only the outer 0.004 Mq is not resolved. Therefore, 
we can say with certainty that this stellar collision product 
does not have a surface convection zone larger than 0.004 
Mq. Onr evolutionary models show that the stellar colli¬ 
sion products do not develop surface convection zones dur¬ 
ing their thermal contraction to the main sequence. 
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Figure 4. Structure of off-axis collision products. The line styles 
are the same as in Figure ^ Again, we see a convergence to the 
fOO 000 particle simulation. 


As well as studying the head-on collisions, we performed 
one off-axis collision with a pericentre separation of 0.25 Rq. 
Figure ^ shows the structure profiles of the collision prod¬ 
ucts at various resolutions. The line styles are the same as 
in Figure except that we have not run a million-particle 
off-axis simulation. In addition to profiles of pressure, tem¬ 
perature, density and hydrogen mass fraction, we show the 
angular velocity and radius profiles. Again, there is clearly 
a convergence of the structural properties as the number of 
particles increases. The 100 000 particle (short dashed line) 
and 300 000 particle (long dashed line) simulations are al¬ 
most indistinguishable. 

The final question we wish to investigate with these high 
resolution simulations involves the material which is thrown 
off during the collision. Between 1 and 10% of the original 
stellar material becomes unbound from the system as it is 
shock-heated during these kinds of collisions, depending on 
t he masses of the par ent stars and the impact parameter 
(Lombardi et al. 1996). Other material is thrown to large 
distances but remains bound and falls back onto the colli¬ 
sion product. If some of this material has sufficient angular 
momentum, it will form a disk around the collision prod¬ 
uct. The existence of such a disk could prove important for 
transporting angular momentum from the central collision 
product. 

Figure ^ is a snapshot of the 300 000 particle off-axis 
stellar collision taken 50 dynamical times after the collision 
began. The star has reached hydrostatic equilibrium and is 
settled down to its post-collision state. The contours show 
the density of the material in a slice through the x-z plane 
of the collision product. The collisions were constrained to 
occur in the x-y plane, so the product is rotating around the 
z axis. It is rotating rapidly, and therefore is highly flattened 


Figure 5. Density profile of an off-axis stellar collision product 
between two 0.6 Mq stars, modelled using 300 000 particles. This 
is a 5x5 Rq box in the x-z plane around the collision product. 
The logarithmically spaced density contours are 0.5 dex apart. 
The collision product is rotating around the z axis, and shows 
significant flattening due to rotation. However, there is no indi¬ 
cation of a disk. 


by rotation. However, there is no indication that a circum- 
stellar disk is present. This is in agreement with the results 
of Benz & Hills (1987), who find that disks are present only 
for simulations with larger impact parameters. 


4 SUMMARY AND DISCUSSION 

We have performed the highest resolution SPH simulations 
of collisions between main-sequence stars to date, using up 
to 10® SPH particles. The collisions were between equal mass 
main-sequence stars in globular clusters - the stars have 
masses of 0.6 Mq, metallicity Z=0.001, age of 15 Gyr and 
a relative velocity at infinity of 10 km/s. We varied the res¬ 
olution of our simulations from 10"^ to 10® SPH particles, 
and we looked at a head-on collision and an off-axis colli¬ 
sion. These simulations took between ~ 6 CPU hours and 
~ 5000 CPU hours, and the CPU time scaled approximately 
linearly with N log N. We had three goals for performing this 
research. 

Our first goal was to determine the optimal resolution 
for SPH simulations of this kind. Since we want to use the 
results of the SPH simulations to study the detailed struc¬ 
ture and evolution of the collision products, we need to have 
an accurate model of the distribution of the SPH particles 
and their properties. On the other hand, we do not want to 
waste valuable computer resources going to a higher reso¬ 
lution than is necessary. We have shown that both head-on 
and off-axis simulations with ~ 100 000 particles give essen¬ 
tially the same information as runs with higher resolution, 
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and therefore we suggest that most simulations with similar 
requirements can confidently use on order 10® particles. 

Our second goal was to use these high resolution simula¬ 
tions to study the outer few percent of the collision products. 
In particular, we were interested in the existence of a surface 
convection zone, for two reasons. One consequence of a con¬ 
vection zone is mixing of elements in the convection zone. If 
the zone reaches deep enough in a star to dredge up nuclear 
processed material, for example, the surface abundances of 
elements like C, N and O will be different from what we 
would expect for primordial material. In blue stragglers and 
other collision products, surface abundances could be used 
to determine, perhaps, the masses of the parent stars. How¬ 
ever, this can only be done if we understand whether the 
abundances should be mixed by convection or not. The sec¬ 
ond consequence of a surface convection zone is angular mo¬ 
mentum loss. S tars can lose a ngular momentum through a 
magnetic wind (Kawaler 19881). This process is most effective 


all off-axis collisions have too much angular momentum), 
or collision products have some way of losing their angu¬ 
lar momentum. The two standard ways stars are thought 
to lose angular momentum is through a magnetic wind, or 
through disk-locking. We have shown that both these sce¬ 
narios are not viable, since these collision products do not 
have surface convection zones or circumstellar disks imme¬ 
diately after the collision. Howev er, it is interesting to note 
the work of Durisen et al. ( 1986|), who show that if a poly¬ 


in stars with deep convection zones, since they can support 
stronger magnetic fields and since the angular momentum 
is drained out of the entire (uniformly rotating) convection 
zone. 


trope is rotating very rapidly, some outer material is thrown 
off into a disk, leaving a stable (but non-axisymmetric) cen¬ 
tral object. In order for this to happen, the star has to have 
a ratio of rotational kinetic energy to gravitational energy 
P > 0.3. At the end of the collision, our off-axis collision 
product has /3 0.045, and is therefore stable against any 

bar instability. However, as the star evolves and contracts, 
it is possible that /3 could becomes high enough to trigger 
some of these hydrodynamic rotational instabilities. 

In this paper, we have limited ourselves to two very spe¬ 
cific collisions. The two stars involved in the collision were of 
equal mass, and we only investigated two choices of impact 
parameter. There have been other studies of main-sequence 


collision investigated in this oaoer does not have a surface 

([Benz & Hills 1987; 

Lombardi et al. 1996 

; Sandquist et al. 

convect 

ion zone larger than 0.004 Mq. For comparison, the 

1997). These studies have clearly and carefully outlined the 


Sun is considered to ha ve a fairly shallow co nvection zone, 
with a mass of 0.02 Mq ([Guenther et al. 1992). A convection 


zone less than about 0.01 Mq is certainly not deep enough 
to substantially modify the surface abundances of any ele¬ 
ment, with the possible exception of lithium, beryllium and 
boron (which are destroyed at low temperatures). Angular 
momentum loss is also going to be very small for stars with 
such small convection zones. 

Our final goal was to determine whether the off-axis 
stellar collision products had a circumstellar disk. If a ro¬ 
tating star has a disk with a mass of about 0.01 Mq, and it 
is locked to the disk by a magnetic field, then angular mo- 
mentum will be transported from the star to the disk by the 


field (Konigl 1991). The net effect of the angular momen¬ 


tum transport will be to slow the star’s rotation rate. This 


effect is thought to operate in p re-main-sequence stars (Sills 
Pinsonileault fc Terndrup 2000|, and references therein). Un¬ 


fortunately, our simulations do not show any evide nce for a 
disk. Based on previous work (Benz & Hills 1987), we ex¬ 


pect that blue stragglers which were formed from collisions 
with larger impact parameters will have circumstellar disks. 
However, the low impact parameter collision we modelled in 
this simulation is still rotating too quickly to contract to the 
main sequence without losing some angular momentum. 

Sills et al. (2001) showed that blue stragglers that are 


formed by stellar collisions have an angular momentum 
problem. The collision products are formed with large to¬ 
tal angular momentum, and some of them are rotating near 
their break-up velocity. The collision products are large ob¬ 
jects, and they start to contract to the main sequence. As 
they contract, they spin up since they have no way to lose 
any angular momentum. They begin to reach break-up ve¬ 
locities and become unstable. Even by shedding material, 
they cannot remove enough angular momentum to remain 
bound, and they tear themselves apart. We are forced to one 
of two conclusions; either blue stragglers are not made by 
stellar collisions (since head-on collisions are very rare, and 


properties of collisions under different circumstances. These 
high resolution simulations were meant to answer some spe¬ 
cific questions about the details of stellar collisions and the 
resulting collision products. We can generalize our current 
results to all main-sequence star collisions in globular clus¬ 
ters under the limitations outlined in the previous papers. 
Therefore, we feel justified in claiming that no collision prod¬ 
ucts have surface convection zones, and that collisions do 
not create long-lived circumstellar disks. The search to un¬ 
derstand the angular momentum evolution of blue stragglers 
will have to turn to other avenues, possibly involving a more 
detailed look at the combination of stellar evolution and hy¬ 
drodynamic instabilities after the collision. 
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